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Abstract 
(N 

> Introduction 

in 

£T} ■ A description of quark bound states demands an application of nonperturba- 

tive approaches in QCD, which comes in the strong coupling regime (as ~ 1) 
within the region of the hadronization at large distances (r ~ 1/Aqcd, 

G\ ; Aqcd ~ 200 MeV). 

^ ■ Considering the bound states with the heavy quarks (nig ^> Aq C d), one 

derives some regularities, simplifying the solution of problem. 
Ph! For the hadrons with a single heavy quark, the virtualities of the heavy 

quark are low (~ Aqcd) and one can make the expansion of the heavy 
quark QCD action over the small parameter Agc^/mg. In the leading ap- 
proximation, the effective action possesses the symmetry with respect to the 
substitution of a heavy quark, moving with a velocity v, by any other heavy 
quark, moving with the same velocity v and having an arbitrary orienta- 
tion of its spin |rj. The symmetry allows one to state the scaling law for 
leptonic constants of the heavy mesons, containing a single heavy quark, 
and the universal dependence of the form factors for exclusive semileptonic 
weak transitions between the heavy hadrons such as B — > D^lv, so that the 
universal function has fixed normalization at the zero recoil point. 

For the heavy quarkonioum (QQ'), the nonrelativistic heavy quark mo- 
tion inside the bound state has allowed one to develop the approach of phe- 
nomenological potential models. The QCD-motivated potentials combine 
the linear rise of the quark interaction energy at the large distances and the 
Coulomb-like interaction at the small distances M. Such models reach the 
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high accuracy in the description of the mass spectra of heavy quarkonia (the 
(cc) charmonium and the (bb) bottomonium), Sm ~ 30 MeV, however, their 
accuracy in the description of leptonic constants is very low (~ 25%). 

Another powerful tool for the description of bound states with the heavy 
quarks has become the QCD sum rules ||, combining perturbative calcula- 
tions and an account of contributions by the vacuum expectation values of 
composite operators, i.e. by the quark gluon condensates such as < mqq >, 
< <y s G 2 > and so on. However, making the consideration with the finite 
number of terms in the QCD perturbation theory for the Wilson's coeffi- 
cients and taking into the account only the restricted set of quark-gluon 
condensates, results of the QCD sum rules get unphysical dependence on an 
external parameter of the sum rule scheme (the number of spectral density 
moments n or the Borel transformation parameter a). This dependence es- 
sentially decreases the predictive power of QCD sum rules. An additional 
parameter is also the threshold s t h, discriminating the resonant region from 
the hadronic continuum. Moreover, the weight functions, rapidly dropping 
with the energy increase, define the averaging scheme of QCD sum rules and 
do not allow one to draw some conclusions on the contributions by the higher 
excitations of quarkonium with the given quantum numbers, so that these 
contributions are practically neglected. 

The QCD sum rule scheme, using the data on the quarkonium mass 

spectrum has been recently offered in ref.0. In the scheme, the conditions 

of nonrelativistic quark motion and the small value of ratio Aqcd/i^q have 

allowed one to state the scaling relation for the leptonic constants of S'-wave 

quarkonia (QQ) 

P 

— = const. , (1) 

independently of the heavy quark flavours. The generalization of eq. ([]]) for 
the heavy quarkonium (QQ') has been considered in ref.||, so that 

P f M V 

M W) = C0USt ' (2) 

where // = mQrriQi j (ttiq + rriQ>) is the reduced mass of (QQ') system. 

In the present paper we use the QCD sum rule scheme of refs. @, |5| to de- 



rive the relation for the leptonic constants of nS'-states of heavy quarkonium 

f = ^, (3) 

that does not depend on the heavy quarkonium content. 

In Section 1 we describe the QCD sum rule scheme, using the spectro- 
scopic data, and derive eq.@. In Section 2 we make the phenomenological 
analysis of relation ([|) and show, that it gives a good description of the 
experimental relations for the leptonic constants in the ip- and T-particle 
families. In the Conclusion we discuss the obtained results. 

1 QCD Sum Rules for Heavy Quarkonium 

Let us consider the two-point correlator functions of quark currents 

IV(g 2 ) = iJd A xe iqx <0\TJ fl (x)Jl(0)\0>, (4) 

U P {q 2 ) = i I ' d A xe iqx < 0|TJ 5 (x)4(0)|0 > , (5) 



where 



J/t{x) = Qi(xh^Q 2 (x) , (6) 

h{x) = Qi(x)-f 5 Q 2 (x) , (7) 

(8) 



Qi is the spinor field of the heavy quark with i = c,b. 
Further, write down 

IV = {-g, v + q -f-y v {q 2 ) + ^n 5 (g 2 ) , (9) 

where ITy and lis are the vector and scalar correlator functions, respectively. 
In what follows we will consider the vector and pseudoscalar correlators: 
IL/(g 2 ) and U P (q 2 ). 

Define the leptonic constants fy and fp 

<0|J Al (x)|y(A)> = ie^f v M v e ikx , (10) 

<0|J 5 ^)l^> = ik„fpe ikx 7 (11) 



where 

JBft{x) = QiixhslnQzi 30 ) > ( 12 ) 

so that 

f A/f 2 
<0|J 5 (x)|P>=2 JP p e ikx , (13) 

mi + r«2 

where |V > and |P > are the state vectors of 1~ and 0~ quarkonia, and 
A is the vector quarkonium polarization, k is 4-momentum of the meson, 

k 2 — M 2 

hpy — lVlpy. 

Considering the charmonium (if), if)' ...) and bottomonium (T, T', T" ...), 
one can easily show that the relation between the width of leptonic decay 
V — > e + e _ and fy has the form 

T(V - e*e-) = f e?al jf , (14) 

where e^ is the electric charge of quark i. 

In the region of narrow nonoverlapping resonances, it follows from eqs.(^) 
- (0) that 

±3roI#"V) = J2fvnM v J(q 2 -M Vn ), (15) 

-^n^V) = J2fpnM Pn ( * , 2 % 2 ~M|J. (16) 

n ^ (m 1 +m 2 ) z 

Thus, for the observed spectral function one has 

-3m<?(g 2 ) = ^mlT^V) + pvA? , ti p) , OT 

TV TV 



where p(q 2 , /i 2 ) is the continuum contribution, which is not equal to zero at 
>/i 2 . 
Moreover, the operator product expansion gives 



q 2 > /i 2 . 



a s 



n (QCD) (g2) = Hipert)^ + ^2) < _A Q 2 > +Q(g2) < ^Q.q, > + ^ 

TV 

(18) 
where the perturbative contribution n^ ert - ) (g 2 ) is labeled, and the nonpertur- 
bative one is expressed in the form of sum of quark-gluon condensates with 
the Wilson's coefficients, which can be calculated in the QCD perturbative 
theory. 



In eq.(P^) we have been restricted by the contribution of vacuum expec- 
tation values for the operators with dimension d = 4. For Cq (q 2 ) one has, 
for instance, || 

{P) = 1 (f_ / 3(3v 2 + l)(l-v 2 ) 2 1_+^ _ 9v 4 + 4v 2 + 3 \ 

G ~ 192m 1 m 2 q 2 \ 2v 5 ^ l-v v 4 ) ' [ ' 

where 

q 2 = q 2 - (mi - m 2 ) 2 , v 2 = 1 - -^^- . (20) 

The analogous formulae for other Wilson's coefficients can be found in Ref. [[|. 
In what follows it will be clear that the explicit form of coefficients has no 
significant meaning for the present consideration. 

In the leading order of QCD perturbation theory it has been found for 
the imaginary part of correlator that || 

QmHy (q 2 ) = - — 2-(3ss — s 2 + 6mim 2 s — 2m\s)6(s — {m\ + rri2) 2 ()2l) 

^mlif rt \q 2 ) = i^( s -( mi - m2 ) 2 )#( s -( mi+m2 ) 2 ), (22) 

where s = s — m\ + m|, s 2 = s 2 — Am^s. 

The one- loop contribution into Qmll(q 2 ) can be included into the consid- 
eration (see, for example, Ref.||). However, we note that the more essential 
correction is that of summing a set over the powers of (a s /v), where v is 
defined in eq . (p0|) and is a relative quark velocity, and as is the QCD interac- 
tion constant. In Ref.|| it has been shown that account of the Coulomb-like 
gluonic interaction between the quarks leads to the factor 

F{V) = IV l-exp(-^) ' (23) 

so that the expansion of the F(v) over as/v <C 1 restores, precisely, the 
one- loop 0(— ) correction 

W-l-^-... (24) 

3 v 

In accordance with the dispersion relation one has the QCD sum rules, which 

state that, in average, it is true that, at least, at q 2 < 

i f^^y ds = 

n J s — q 2 



where the necessary subtractions are omitted. ^hnli^ ha ^(q 2 ) and H^ CD \q 2 ) 
are defined by eqs.(^) - (|P7D and eqs.(^) - fl24|) , respectively, eg . (p5|) is the 



base to develop the sum rule approach in the forms of the correlator function 
moments and of the Borel transform analysis (see Ref.||). The truncation 
of the set in the right hand side of eq.(p5|) leads to the mentioned unphysical 
dependence of the fpy values on the external parameter of the sum rule 
scheme. 

Further, let us use the conditions, simplifying the consideration due to 
the heavy quarkonium. 

1.1 Nonperturbative Contribution 

We assume that, in the limit of the very heavy quark mass, the power cor- 
rections of nonperturbative contribution are small. From eq. (|i~9"D one can see 
that, for example, 

4 P V) « 0(-}—) , A/m lj2 «l, (26) 

where v is fixed, q 2 ~ (m\ + rn.2) 2 , when 9mll' pert ' (q 2 ) ~ (mi + mq) 1 . It is 
evident that, due to the purely dimensional consideration, one can believe 
that the Wilson's coefficients tend to zero as l/m\ 2 . 

Thus, the limit of very large heavy quark mass implies that one can 
neglect the quark-gluon condensate contribution. 

1.2 Nonrelativistic Quark Motion 

The nonrelativistic quark motion implies that, in the resonant region, one 
has, in accordance with eq . (p0|) , 



v -»• . (27) 

So, one can easily find that in the leading order 



Zmlir rt, (s) « %mll ( r rt \s) - ^- 2 s (j£j , (28) 



so that with account of the Coulomb factor 



F(v) ~ ^ ^ , (29) 
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one obtaines 



^ rt) (^| S g) 2 . (30) 



1.3 "Smooth Average Value" Scheme of the Sum Rules 

As for the hadronic part of the correlator, one can write down for the narrow 
resonance contribution 



niTV) = /^E«*-m^), (si) 

J O U ,v> 



nrv) = i— 2 Y,fin ( M r n , 2 ^-mij, (32) 

J s - Q n ( m l + m 2) 

The integrals in eqs.(|3T|)-(|3"2"|) are simply calculated, and this procedure is 
generally used. 

In the presented scheme, let us introduce the function of state number 
n(s), so that 

n{ml) = k . (33) 

This definition seems to be reasonable in the resonant region. Then one has, 
for example, that 



-smU^ s \s) = sf Vn{s) Ts T,0(s- Ml k ) . (34) 

Further, it is evident that 

sp(-«J) = ^ !?»("-*)■ (35) 



and eq.fl31"D can be rewritten as 

n ^V) = /^<„ w ^ ?»(»-*)■ < 36 > 



The "smooth average value" scheme means that 

<" ) tf)=<i^-«)>J 7 ^ r fv^- (37) 
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It is evident that, in average, the first derivative of step-like function in the 
resonant region is equal to 

<-^-£0(n-Jfe)>~l. (38) 

k 

Thus, in the scheme one has 

<IirV)> - j^fvni/-^ 1 , (39) 

<IirV)> - / ds s2f ^(s) dn^_ (40) 

J s — q 2 (mi + m 2 y as 

Eqs.(^S])-(^) give the average correlators for the vector and pseudoscalar 
mesons, therefore, due to eg. (pop we state that 



3m < n^V) >= 3mn( QCD V) , (41) 

that gives with account of eqs.(pO[), (|39"D and fl40|) at the physical points 

Jl = as_ dNU /4/i\ 2 

M n " 7i dn \Mj ' l ' 

where in the limit of heavy quarks we use, that for the resonances one has 

m 1 + m 2 ~ M , (43) 

so that 

f Vn ~ f Pn = f n . (44) 

Thus, one can conclude that for the heavy quarkonia the QCD sum rules 
give the identity of fp and fy values for the pseudoscalar and vector states. 

Eq.(^) differs from the ordinary sum rule scheme because it does not 
contain the parameters, which are external to QCD. The quantity dM n /dn 
is purely phenomenological. It defines the average mass difference between 
the nearest levels with the identical quantum numbers. 

Further, as it has been shown in ref.||, in the region of average distances 
between the heavy quarks in the charmonium and the bottomonium, 

0.1 fm < r < 1 fm , (45) 



the QCD-motivated potentials allow the approximation in the form of loga- 
rithmic law J7| with the simple scaling properties, so 

dn 

const. , (46) 



dM„ 



i.e. the density of heavy quarkonium states with the given quantum numbers 
do not depend on the heavy quark flavours. 

In ref . H it has been shown, that relation ( f4"6D is also practically valid for 
the heavy quark potential approximation by the power law (Martin potential) 
[[J, where, neglecting a low value of the binding energy for the quarks inside 
the quarkonium, one can again get eg. (|46|) . 

In ref. [[D it has been found, that relation (|46|) is valid with the accuracy 
up to small logarithmic corrections over the reduced mass of quarkonium, if 
one makes the quantization of S"-wave states for the quarkonium with the 
Martin potential by the Bohr-Sommerfeld procedure. 

Moreover, with the accuracy up to the logarithmic corrections, as is the 
constant value. Thus, as it has been shown in refs.^j, H], for the leptonic 
constants of S-w&ve quarkonia, the scaling relation takes place 

f 2 fM\ 2 

M U = m " St • < 47) 

independently of the heavy quark flavours. 

Taking into the account eqs.(f43D and fl44|) and integrating eqs. (|39"D , (|40|) 
by parts, one can get with the accuracy up to border terms, that one has 

df n dn a s f^\ 2 

- 2fn ^dM n n = V Mn {M:) ■ (48) 



(49) 



Comparing eqs.p^) and (p8|), one finds 

dfn = 1 

f n dn 2n 

that gives, after the integration, eq.(f|): 

/ni _ n 2 
fn 2 n l ' 
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Table 1: The experimental values of leptonic constants (in MeV) for the 
nS'-bottomonia in comparison with the estimates of present model. 



quantity 


exp. 


present 


h 


715 ± 15 


input 


h 


487 ±16 


506 ± 10 


h 


429 ± 14 


412 ±8 


h 


320 ± 30 


358 ±7 


h 


369 ± 46 


320 ±7 


h 


240 ± 30 


292 ±6 



Relation (^) leads to that the border terms, which have been neglected in 
the writing of eg. (|48|) , are identically equal to zero. 

Thus, under the conditions of small Aqcd/^^q ratio value and the nonrel- 
ativistic heavy quark motion, the described QCD sum rule scheme with the 
"smooth average value" allows one to derive the scaling expression, relating 
the leptonic constants of different nS-wave levels, independently of the quark 
content of heavy quarkonium. 



2 Analysis of Scaling Relation 



First, note that eq . (f47|) , relating the leptonic constants of different quarkonia, 
turns out to be certainly valid for the quarkonia with the hidden flavour (cc, 
bb), where 4/z/M = 1 (see [Q). In that case, to estimate the constant value 
in the right hand side of eq.fl4"?|) we have supposed in refs.j|, |[, that as has 
the value, defining the Coulomb part of potential in the Cornell model (see 
[0]), and the < dM/dn > value is equal to the average distance between the 
nearest S'-wave levels in the bottomonium. 

Second, eq. (^7J) gives estimates of the leptonic constants for the heavy B 
and D mesons, so these estimates are in a good agreement with the values, 
obtained in the framework of other schemes of the QCD sum rules || . 

These two facts show that the offered scheme can be quite reliably applied 
to the systems with the heavy quarks. 

Taking a value of the 15*- level leptonic constant as the input one, we 
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calculate the leptonic constants of higher nS'-excitations in the charmonium 
and the bottomonium. 

The results are presented in Tables [[J |2| and on Figures [I], |2|. One can see 
that eq.(|3]) is in a good agreement with the experimental values of leptonic 
constants for the nS-wave levels of heavy quarkonia || . 

One has to note, that for the fy^s) value we have taken 



if 



i/>(3S) 



/V(3770) + it/. (4040) 



(50) 



since, as it generally accepted, the ip(3770) and ^(4040) states are the results 
of 3D- and SS'-levels mixing in the charmonium. 

Conclusion 



In the framework of the QCD sum rules for the leptonic constants of the heavy 
quarkonia one uses the conditions of low Aqcd/wq ratio value and the non- 
relativistic quark motion in the phenomenological potential, possessing the 
simple scaling properties, one takes into the account the Coulomb-like as/v 
corrections, and in the scheme of the "smooth average value" one derives the 
scaling expression, relating the leptonic constants of rzS'-wave quarkonium 
levels, so 

fn 2 n ± ' 

independently of the heavy quark flavours. 



Table 2: The experimental values of leptonic constants (in MeV) for the 
nS'-charmonia in comparison with the estimates of present model. 



quantity 


exp. 


present 


h 


410 ± 14 


input 


h 


283 ± 14 


290 ± 10 


h 


205 ± 20 


237 ±8 


h 


180 ± 30 


205 ±7 


h 


145 ± 15 


183 ±6 
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Figure 1: The calculated dependence oinS - bottomonium leptonic constants 
and the experimental values of fr( n S)- 

The obtained relation is in a good agreement with the experimental val- 
ues of leptonic constants for the charmonium and the bottomonium, and it 
reflects a small variation of the heavy quark kinetic energy with respect to 
the heavy quark flavours. 
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